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ABSTRACT 

Background: Obesity is a complex trait with both 
environmental and genetic contributors. Genome-wide 
association studies have identified several variants that 
are robustly associated with obesity and body mass 
index (BMI), many of which are found within genes 
involved in appetite regulation. Currently, genetic 
association data for obesity are lacking in Africans— a 
single genome-wide association study and a few 
replication studies have been published in West Africa, 
but none have been performed in a South African 
population. 

Objective: To assess the association of candidate loci 
with BMI in black South Africans. The authors focused 
on single nucleotide polymorphisms (SNPs) in the 
FTO, LEP, LEPR, MC4R, NPY2R and POMC genes. 
Design: A genetic association study. 
Participants: 990 randomly selected individuals from 
the larger Birth to Twenty cohort (a longitudinal birth 
cohort study of health and development in Africans). 
Measures: The authors genotyped 44 SNPs within the 
six candidate genes that included known BMI- 
associated SNPs and tagSNPs based on linkage 
disequilibrium in an African population for FTO, LEP 
and NPY2R. To assess population substructure, the 
authors included 18 ancestry informative markers. 
Weight, height, sex, sex-specific pubertal stage and 
exact age collected during adolescence (13 years) 
were used to identify loci that predispose to obesity 
early in life. 

Results: Sex, sex-specific pubertal stage and exact 
age together explain 14.3% of the variation in log(BMI) 
at age 13. After adjustment for these factors, four 
SNPs were individually significantly associated with 
BMI: FTO rs1 781 7449 (p=0.022), LEP rs1 09541 74 
(p=0.0004), LEP rs6966536 (p=0.012) and MC4R 
rs1 778231 3 (p=0.045). Together the four SNPs 
account for 2.1 % of the variation in log(BMI). Each risk 
allele was associated with an estimated average 
increase of 2.5% in BMI. 

Conclusions: The study highlighted SNPs in FTO and 
MC4R as potential genetic markers of obesity risk in 
South Africans. The association with two SNPs in the 
3' untranslated region of the LEP gene is novel. 



ARTICLE SUMMARY 



Article focus 

■ This is a replication study aiming to reproduce 
BMI association findings from European cohorts 
in a South African population. 

■ This study focused on genes linked to appetite 
control that were previously reported to show 
association with BMI or obesity and included 
FTO, LEP, LEPR, MC4R, NPY2R and POMC. 

■ Adolescent data were used to facilitate the 
identification of genetic loci that predispose to 
obesity early in life, as it is known that 
overweight/obese children have an elevated risk 
of becoming obese adults. 

Key messages 

■ We found four SNPs were individually signifi- 
cantly associated with BMI: FTO rs1 781 7449 
(p=0.022), LEP rs1 09541 74 (p=0.0004), LEP 
rs6966536 (p=0.012) and MC4R rs1 778231 3 
(p=0.045). 

■ Together the four SNPs account for 2.1% of the 
variation in log(BMI). 

■ We also demonstrated that an accumulation of 
risk alleles is linked to a significant increase in 
BMI— individuals with seven risk alleles had an 
11.0% increase in median BMI compared with 
those with two risk alleles. 

Strengths and limitations of this study 

■ This study provides the first preliminary evidence 
of the role of genetic variants in obesity risk in an 
adolescent black South African population. 

■ This study was only moderately powered to 
detect association with BMI, and not all genes 
were exhaustively investigated. 

■ TagSNP selection would have been enhanced if 
South African data were available for this 
approach. 
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Candidate genes associated with BMI in South Africans 



INTRODUCTION 

Being overweight or obese is a significant risk factor for 
the development of chronic diseases like type II diabetes 
and cardiovascular disease and is becoming increasingly 
common in low- and middle-income countries. 1 Both 
environmental and genetic factors influence body mass, 
with the heritability of body mass index (BMI, kg/m 2 ) 
estimated at 40%-70%. 2 

The study of syndromic obesity and genome-wide asso- 
ciation studies (GWAS) of common measures of adiposity 
(eg, BMI) have provided evidence for genetic risk loci for 
obesity. 3 4 Current GWAS have primarily focused on 
populations of European origin, and only one GWAS for 
anthropomorphic traits, including BMI, has been 
performed in an indigenous African population. This 
study was carried out in a cohort of 1188 Nigerian subjects 
and replicated the association of MCAR with BMI, but no 
other significant associations were detected. 5 However, the 
study was underpowered to detect modest effects. 
Although several replication studies have been performed 
in African— Americans, 6 7 only a few are reported in well- 
characterised African populations 8 9 and showed some 
supporting evidence of the role of FTO in a Nigerian 
population but not in the Gambia. Genetic association 
studies in African populations have the significant advan- 
tage that linkage disequilibrium (LD) generally exists over 
a shorter genomic distance, potentially increasing the 
efficiency of the identification of causal variants. 10 

Several of the likely causal genes for increased body fat 
mass are involved in appetite regulation. Specifically, 
constituents of the leptin— melanocortin neurotransmis- 
sion pathway have been implicated in both monogenic 
and polygenic obesity. Recent evidence that FTO (the only 
gene consistently associated with measures of obesity) acts 
on the central nervous system-mediated control of food 
intake 11 lends further support for the hypothesis that 
obesity may be a heritable neurobehavioral disorder. 

Obesity levels are rising at an alarming rate among 
black South Africans, with 29% and 57% of men and 
women, respectively, being overweight or obese. The 
major impact of obesity on the development of chronic 
diseases within South Africa is demonstrated by data 
showing that in the year 2000, 87% of type II diabetes 
cases and 68% of hypertensive disease were attributable to 
a BMI of &21. 1 The multifactorial nature of these 
conditions creates a challenge, and understanding the 
genetic contribution to obesity in black South Africans 
may influence the nature of effective interventions. The 
aim of this study was therefore to assess the association of 
genetic variants with BMI in black South African adoles- 
cents. We focused on genes linked to appetite control that 
were previously reported to show association with BMI or 
obesity: FTO, LEP, LEPR, MC4R, NPY2R and POMC. 

METHODS 
Participants 

The investigation focused on participants of the South 
African Birth to Twenty (Bt20) cohort study. Bt20 is 



a longitudinal study of the health and well-being of 
children who were born in the metropolitan area of 
Soweto, Johannesburg, during a 6-week period in 1990. 
Participants in Bt20 were enrolled at birth, and detailed 
information has been collected from the participants 
and their care givers over the past 20 years. 12 Following 
informed consent, blood samples were collected and 
DNA extracted using a salting-out method. 13 The Bt20 
research programme has received clearance by 
the Human Research Ethics Committee: (Medical) of 
the University of the Witwatersrand (M010556). The 
Federal-Wide Assurance registration number of the 
committee is FWA00000715. 

A subset of individuals from the Bt20 cohort (524 
female and 466 male adolescents, total n=990) were 
randomly selected for this study. This represents 43% of 
the total cohort (2290 participants). The following 
calculation illustrates that this study was adequately 
powered. Group sizes of 748, 204 and 21 (total=973), 
corresponding to rsl0954174 genotype frequencies 
observed in this study, will achieve 93% power to detect 
differences among the means of BMI using an F test with 
a 0.05 significance level, as is standard in linear regres- 
sion models. The size of the variation in the means (say 
18, 19 and 20) is represented by their SD, which is 0.47. 
The common SD within a group is assumed to be 4. 

Adolescent data (participants were 13 years of age at the 
time measurements were done) were used to facilitate the 
identification of genetic loci that predispose to obesity 
early in life, as it is known that overweight/ obese children 
have an elevated risk of becoming obese adults. 14 Anth- 
ropometric measurements were obtained using standard 
methods. 13 BMI was computed as weight (measured in 
kilogram) divided by the square of the height (measured 
in metres) of an individual. Pubertal stage was assessed 
using a validated self-assessment method drawing upon 
the Tanner Sexual Development Tool, 16 and exact age at 
time of measurement was recorded. 

Population substructure 

A set of 18 ancestry informative markers (AIMs) were 
included to assess population substructure within the 
study group, based on previously published data and 
local unpublished data. 18 Population substructure was 
assessed computationally using STRUCTURE V.2.3.3 
software. 19 For comparative purposes, genotypes of the 
AIMs were obtained from HapMap (release 24) for three 
additional populations — the Yoruba from Nigeria (YRI), 
the Han Chinese population (CHB) and Utah residents 
with European ancestry (CEU). STRUCTURE imple- 
ments the value K to distinguish parental populations, 
with K referring to the number of inferred population 
clusters. The true value of K may not always be known, 
but the smallest value of K that captures the majority of 
structure in the data should be obtained. Ten iterations at 
K=l to K=4 with a burn-in of 30 000 followed by 100 000 
iterations were implemented. The K value with the 
highest average posterior likelihood score was calculated 
and used to identify the best cluster assignment. 



2 



Lombard Z, Crowther NJ, van der Merwe L, etal. BMJ Open 201 2;2:e000873. doi:1 0.1 1 36/bmjopen-201 2-000873 



Candidate genes associated with BMI in South Africans 



Single nucleotide polymorphism selection 

Single nucleotide polymorphisms (SNPs) previously 
shown to be associated with BMI or obesity in FTO, LEP, 
LEPR, MC4R, NPY2R and POMC were identified from 
the literature. In addition, tagSNPs were selected for 
LEP, NPY2R and POMC to ensure that known common 
polymorphisms are either directly assayed or associated 
with a tagSNP in these genes. SNPs from the literature 
were not force included in tagSNP selection. The 
Tagger algorithm 20 was used to select tagSNPs in 
a multimarker approach at r 2 >0.80 and with a minor 
allele frequency >0.05 among publicly available African 
data (YRI) from the HapMap data set (release 24). The 
Illumina Assay Design Tool (Illumina, San Diego, Cali- 
fornia, USA) was used to assess the SNPs' compatibility 
for inclusion in a custom-design GoldenGate VeraCode 
Assay, and 62 SNPs were selected for genotyping 
(supplementary table 1). 

Genotyping 

Genotyping was performed using the Illumina Golden- 
Gate VeraCode assay (Illumina). Quality control was 
performed on all raw genotype data according to the 
supplier's specifications using BeadStudio V.2.0 software 
(Illumina) (summarised in Tindall et af 1 ). Standard 
quality control filters were applied to the data 22 : minor 
allele frequency >0.01, SNP missingness rate <0.05, 
individual missingness rate <0.2 and Hardy— Weinberg 
equilibrium (HWE) <1X10" 4 . 

LD plots 

LD was visualised using Haploview V.4.2, 23 implementing 
the CI method of Gabriel et af 4 to construct haplotype 
blocks. The chromosomal location of SNPs was obtained 
from NCBI (genome build 37.1). 

Statistical analysis 

The programming environment R, and R package 
genetics were used for all statistics (http://www.r-project. 
org). 

Individual genotype associations with BMI 

Genotype and allele frequencies were calculated 
through direct counting. The distribution of BMI 
was skewed to the right, and therefore, BMI was log- 
transformed to approximate normality for all analyses. 
Linear models were used for all analyses (individual 
association, joint association, sex-specific association and 
risk score association), starting by comparing BMI 
between genotype groups, enabling adjustment for 
confounders. As BMI correlated significantly with sex, 
sex-specific pubertal stage and exact age, some analyses 
were adjusted for these variables by including them in 
the linear models as covariates. All association p values 
are from linear models of log(BMI). Exact tests of HWE 
were also performed. Because of its skewed distribution 
in this study group, BMI is summarised as median values 
and IQRs, with lower quartile and upper quartiles 
reported. 



Furthermore, each genotype was modelled additively 
as the number of minor alleles present. Because we 
modelled log (BMI), the modelled effect of each allele is 
reported as a percentage change (increase or decrease) 
in BMI. 

Joint model and allele risk score of all associated SNPs 

We also created and summarised a linear model 
containing all the variants that were individually signifi- 
cantly associated with BMI. We used all significantly 
associated variants to calculate a risk allele score by 
counting the number of alleles associated with 
increasing BMI in each adolescent. This risk allele score 
was also assessed with a linear model, independent of the 
three confounders. 

Sex-specific effects 

In order to detect sex-specific effects, we modelled the 
interaction between sex and each variant and summar- 
ised the significant results. 

Correcting for multiple testing 

Correcting for multiple testing is a contentious issue, 
and some approaches (such as Bonferroni correction) 
are considered overconservative and one risks the 
rejection of true findings. 25 26 Given the strong prior 
information about the role of the variation tested here in 
obesity, we considered this a replication study, and 
therefore, p values below 0.05 were considered signifi- 
cant. For tagSNPs, tests of associations could be consid- 
ered discovery rather than replication, but since these 
markers are correlated due to LD, the Bonferroni 
assumption of independence is not upheld. 

RESULTS 
Participants 

The study group consists of 524 (53.0%) female and 466 
(47.0%) male adolescents, with a mean age of 13.7 years 
(SD=0.2 years). Female adolescents had a higher 
median BMI (19.4 (IQR=17.8-22.3)) than male 
adolescents (18.0 (IQR=16.6-19.7); p<0.0001). BMI 
increased across the five pubertal stages in both male 
(p<0.0324) and female adolescents (p<0.0001), but at 
significantly different rates. Exact age (days) correlated 
positively with BMI in male adolescents (p=0.027) but 
not in female adolescents (p=0.305). In a joint model, 
sex, sex-specific pubertal stage and age each correlated 
independently with BMI and together explained 14.3% 
(p<0.001) of the variation in log(BMI). Pubertal effect 
when assessed alone had the greatest impact (10.7%) 
followed by sex (5.7%). As a result, all analyses were 
adjusted for age, sex and pubertal stage per sex. Unad- 
justed p values of individual tests are presented for 
comparison (table 1). 

Population substructure 

One AIM failed quality control measures and was 
subsequently excluded from analysis. Analysis of allele 
frequencies of AIMs in the South African cohort 
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(supplementary table 2), and comparison to three other 
ethnic groups (supplementary table 3), confirmed that 
there was no significant population substructure within 
the study cohort (supplementary figure 1). The South 
African population and the YRI group showed distinctive 
clustering away from the Caucasian and Asian popula- 
tions. Although the South African cohort and the YRI 
had similar clustering patterns, higher order clustering 
of all individuals included in this analysis (ie, higher 
values of K, where K indicates the number of parental 
populations) was shown to highlight the genetic 
distinction between these two African groups. 

LD plot 

To illustrate the unique patterns of LD in the South 
African population, an LD plot of the gene most signif- 
icantly associated with BMI (LEP) was constructed 
(supplementary figure 2). Evidence of three LD blocks 
covering the gene is observed, with the 3 end of the 
gene in particular exhibiting high LD. 

Statistical analysis 

Individual genotype associations with BMI 

The genotype and allele frequencies of all candidate 
SNPs, as well as p values for tests of additive allelic 
association with log(BMI) are shown in table 1. All SNPs 
were in HWE (table 1). Table 2 presents the median and 
IQR of BMI for SNPs that were significantly associated 
with BMI following adjustment for sex, pubertal stage 
and age. Estimated effect sizes (percentage difference in 
BMI associated with each additional minor allele 
compared with the major allele homozygote) adjusted 
for age, sex and pubertal stage are also presented. For 
two SNPs, LEP rs6966536 and FTO rsl 781 7449, the 
minor allele (G in both cases) is associated with 
increased BMI, whereas the major allele is associated 
with increased BMI for the other two SNPs, LEP 
rsl0954l74 and MC4R rsl7782313 (G and T, respec- 
tively). When referring to 'risk' allele, it is the allele that 
is associated with increased BMI. 

For LEP rs6966536, each G-allele was associated with 
a 2.9% higher BMI (p=0.012), whereas for LEP 
rsl 09541 74, each A-allele is associated with an estimated 
decrease of 4.0% in BMI (p<0.001). Thus, for an 
adolescent of the same age, sex and pubertal stage, we 
would expect an individual with the A/G genotype at 



rsl0954174 to have a BMI that is 4.0% less than that of 
a G/ G homozygote; and the A/A homozygote is 
expected to have an 8.0% lower BMI compared with the 
G/G homozygote. Adolescents who are A/G heterozy- 
gotes for LEP rs6966536 are expected to have a 2.9% 
higher BMI than A/A homozygotes, while the G/G 
homozygotes would have a 5.8% higher BMI. 

Sex-specific effects 

Two SNPs — MC4R rsl2970134 and NPY2R 
rsl 1099992 — showed significant sex differences in their 
effects on BMI (results not shown). In female adoles- 
cents, each G-allele of the NPY2R rsl 1099992 SNP 
predicted a decrease in BMI by 3.7% (p=0.003), while in 
male adolescents, the effect was not significant 
(p=0.523). In male adolescents, each A-allele of the 
MC4R rsl2970134 SNP predicted an increase in BMI of 
5.2% (p=0.011), while in female adolescents, the effect 
was not significant (p=0.785). These SNPs were not 
associated with log(BMI) in the group as a whole 
(table 1). 

Joint model 

To assess the combined effect of the four SNPs that were 
individually associated with log(BMI) and the three 
confounders (sex, pubertal stage and age), a model was 
built based on 908 participants who had complete 
information. The model explained 16.4% of the varia- 
tion in log(BMI), so that the four variants contributed 
2.1% to the variation, independent of the three 
confounders (table 3) .The two L£PSNPs had the largest 
effects, followed by the FTO and then the MC4R SNPs, 
which each had similar effect sizes (2.2 and —2.0, 
respectively). The seven FTO SNPs in this study 
explain 0.6% of variation in log(BMI) after adjusting for 
age-, gender- and sex-specific pubertal stage and 1.4% 
unadjusted. 

Allelic risk score 

Figure 1 shows the number of adolescents with each 
number of additive risk alleles (the allele associated with 
higher BMI — four loci and therefore a maximum of 
eight risk alleles in an individual) as well as the median 
and IQR of BMI for each group. Since there were <10 
individuals each in the categories of one risk allele 
(seven individuals) and eight risk alleles (one 



Table 2 Summary BMI statistics for SNPs that were individually associated with BMI 















Median (IQR) of BMI per genotype 








Gene 


SNP 


Typed 


Modelled 


MA 


MAF 


mm* 


mM* 


MM* 




Effect 


FTO 


rsl 781 7449 


977 


952 


G 


0.37 


18.5 (16.9-20.7) 


18.7 (17.2-21.1) 


19.1 (17.6- 


-22.5) 


1.9 


LEP 


rsl 09541 74 


973 


948 


A 


0.13 


18.8 (17.1-21.3) 


18.4 (17.2-20.5) 


17.4 (16.8- 


-19.1) 


-4.0 


LEP 


rs6966536 


972 


948 


G 


0.14 


18.6 (17.1-20.8) 


18.9 (17.1-21.5) 


19.2 (18.5- 


-21.6) 


2.9 


MC4R 


rsl 778231 3 


969 


945 


C 


0.25 


18.9 (17.4-21.2) 


18.5 (17.0-20.9) 


19.0 (16.5- 


-20.6) 


-1.8 



Effect sizes are from individual models, estimated individual percentage change in BMI for each minor allele, independent of age, sex and sex- 
specific pubertal stage. Effect sizes correspond to p values (adjusted) in table 1 . 
*m denotes major allele; M is minor allele. 

BMI, body mass index; MA, minor allele; MAF, minor allele frequency; SNP, single nucleotide polymorphism. 
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Table 3 Joint model of SNPs individually associated with 
log(BMI) 



Gene 


SNP 


MA 


Effect 


p Value 


FTO 


rs1 781 7449 


G 


2.2 


0.007 


LEP 


rs1 09541 74 


A 


-3.2 


0.005 


LEP 


rs6966536 


G 


2.7 


0.020 


MC4R 


rs1 778231 3 


C 


-2.0 


0.028 



The model explains 16.4% of the covariation in log(BMI). 
p Values are from the joint model, so they are adjusted for age, sex 
and sex-specific pubertal stage and all other SNPs in the model. 
Effect = estimated % change in BMI of each minor allele, 
independent of other confounding variables (age, sex and sex- 
specific pubertal stage) and all other SNPs in the model (908 
adolescents included in model). 

BMI, body mass index; MA, minor allele; SNP, single nucleotide 
polymorphism. 



individual), these were excluded from the figure. The 
estimated increase in BMI, for each risk allele, is 2.5%, 
independent of age, sex and sex-specific pubertal stage, 
and individuals with seven risk alleles had an 11.0% 
increase ((20.0— 18.0)/18.0) in median BMI compared 
with a group of individuals with two risk alleles. 

DISCUSSION 

This study provides preliminary evidence of the role of 
genetic variants in obesity risk in an adolescent black 
South African group. Individual SNPs in FTO, LEP and 
MC4R genes were shown to be associated with an 
increase in BMI, and the cumulative effect of the risk 
alleles (one to eight alleles) in FTO, LEP and MC4Rwas 
demonstrated. Together they contribute 2.1% to the 
variation observed in log(BMI) in this group of subjects, 
after adjusting for age, sex and sex-specific pubertal 
stage. Recent GWAS in populations of European 
ancestry have identified 32 common loci associated with 
BMI, explaining only 1.5% of the interindividual varia- 
tion in BMI. 4 We also demonstrated that an accumula- 
tion of risk alleles is linked to a significant increase in 
BMI (figure 1). In our study, individuals with seven risk 



alleles had an 11.0% increase in median BMI compared 
with individuals with two risk alleles. 

We provide preliminary evidence that SNPs in the 3 
untranslated region (3'UTR) of LEP (rsl0954l74, 
p=0.0004 and rs6966536, p=0.012) are associated with 
BMI in South Africans. Leptin is an integral molecule in 
energy homeostasis, and circulating leptin is primarily 
produced by adipocytes. The level of circulating leptin 
acts as a satiety signal. Coding or splice site mutations in 
the LEP gene render the leptin protein unable to signal 
through its receptor, which leads to severe childhood- 
onset obesity. 27 28 These mutations are rare in the 
general population, and their impact on the common 
obesity phenotype remains unclear. 

The LEP SNPs that have been most consistently asso- 
ciated with obesity are rs2167270 (also known as LEP 
19G>A) in the untranslated first exon, 29-31 and 
rs7799039 (LEP -2548G>A) in the 5 -region of the 
gene. 32-35 Studies in a Brazilian population 34 36 and 
adult Samoans 37 have also shown obesity to be correlated 
with variants of the LEP 3 hypervariable region (HVR). 
The 3 ETVR is a (TTTC) n tetranucleotide repeat located 
approximately 2000 bp downstream from LEP. The 
3 TJTR SNPs associated with BMI in this study do not 
share a haplotype block with the 3 ETVR in Africans and 
could be a manifestation of a uniquely African associa- 
tion to BMI and obesity. There are several plausible 
mechanisms by which a 3 UTR variant could impact on 
the functioning of a gene. They could reside in regions 
of miRNA binding sites or may affect gene expression in 
other ways as 3 UTRs have been shown to be involved in 
posttranscriptional regulation of gene expression 
through the modulation of nucleocytoplasmic mRNA 
transport, translation efficiency, subcellular localisation 
and message stability. 38 

Variants in FTO (the only gene consistently associated 
with obesity) have been associated with adiposity 
measures in Europeans, 4 39 Asians 40 41 and Afri- 
can—Americans, 7 with the strongest effect seen with SNP 
rs9939609. Our study showed a modest association 



Figure 1 Observed cumulative 
effect of high-risk alleles on body 
mass index (BMI) in a black South 
African adolescent cohort. The 
four risk alleles included in this 
analysis are: LEP rs6966536 (G- 
allele), LEP rs1 09541 74 (G-allele), 
MC4R rs1 778231 3 (G-allele) and 
FTO rs1 781 7449 (T-allele). Bar 
chart (left axis) shows number of 
adolescents with specific number 
of obesity risk alleles. Points show 
corresponding median (right axis) 
and vertical lines represent the 
IQR of BMI. 
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between the FTO SNP, rsl7817449, and BMI after 
adjustment for sex, age and sex-specific pubertal stage. 
An estimated effect size of a 1.9% increase in BMI for 
each minor allele of rs 178 17449 was observed. Inconsis- 
tent associations between FTO variants and obesity 
measures have been observed in native African popula- 
tions. In Gambians, no associations between FTO and 
weight-for-height z scores were observed, rsl 781 7449 was 
directly assayed in that study but showed no association. 
It should be noted that the Gambia has a very low obesity 
rate, which could possibly explain the null result. In West 
Africans, rsl 78 17449 was not associated with BMI, but 
several other FTO variants were found to be associated. 8 
The role of genetic variation at the FTO locus in predis- 
posing to obesity in African populations warrants further 
investigation, specifically in relation to the epidemiolog- 
ical transition and access to a calorie-rich diet. 

The rsl 77823 13 SNP is one of several variants in MC4R 
linked to BMI in European populations and has been 
shown to contribute to a significant increase in extreme 
obesity in adults. 42 In the present study, this variant 
showed a borderline association with BMI (p=0.045), 
suggesting that this SNP is not the functional variant. It 
is located 3 to the gene, and it is therefore possible that 
it may be in LD with another polymorphism that disrupts 
the function of cis-acting elements involved in the 
control of MC4R gene transcription. 43 

In this study, the SNP rsl 1099992 in NPY2R was 
observed to associate with BMI in female subjects only 
(p=0.003), whereas, MC4R rsl 29701 34 showed a sex- 
specific effect in male subjects (p=0.011). It has been 
observed that SNPs in the NPY2R gene have effects on 
obesity in male subjects only. 44 46 One of these studies 
also showed that a SNP in PYY (encoding a ligand of 
NPY2R) associates with obesity-related traits in female 
subjects only. 45 Furthermore, mutations in the MCAR 
gene that have been associated with morbid obesity have 
stronger effects in female subjects than male subjects. 47 
These data suggest that the phenotypic effect of 
sequence variation in the MCAR gene is differentially 
modified by sex depending on the site of the variant. 
Thus, genetic variants relevant to appetite-regulatory 
pathways seem to have sex-specific effects on body 
composition, but the mechanism by which this occurs is 
still poorly understood. 

The data in this study are derived from a cohort of 
adolescents in the midst of puberty. It is therefore 
possible that the effects on weight of some poly- 
morphisms may have been masked by puberty-associated 
changes in body fat mass. 48 0 It is also possible that the 
effects of some polymorphisms on BMI are observed only 
later in life. 51 However, it has been shown that the FTO 
and MCAR polymorphisms do influence anthropometry 
in adults and children and that other polymorphisms 
identified in GWAS in adults also have observable effects 
on childhood measures of adiposity. 42 52-54 

Replication has become the gold standard for 
assessing the validity of statistically significant results 



from association studies. However, a true result can fail 
to replicate in another independent sample for 
numerous reasons including differences in allele 
frequencies and haplotype structure that could influ- 
ence the indirect detection of functional poly- 
morphisms. A recent study showed that the power to 
replicate the statistically significant independent main 
effect of one polymorphism can drop dramatically 
with a change of allele frequency of <0.1 at a second 
interacting polymorphism. Differences in allele 
frequency could also result in a reversal of allelic effects 
where a putative protective allele becomes associated 
with increased risk in a replication study. 55 These 
are particularly important considerations when under- 
taking genomic research and replication studies in 
African populations, with known high levels of genome 
diversity and population substructure between 
different African groups. 56 57 The AIM data obtained in 
this study highlight that there are not high levels 
of heterogeneity within the South African group, but 
that the group differs from other African tribes, such 
as the Yoruba. Despite these differences, 18 AIMs 
cannot resolve the detailed population differences 
between this South African population and other 
African groups, as the Bantu expansion occurred 
rather recently in history (~ 4000 years ago), and one 
would expect that only major differences would be 
observable when interrogating a small number of 
markers. 57 

Further studies on other candidate obesity loci in 
African populations will provide a better understanding 
of the role of variants in these genes in a population with 
a different LD structure to that of Caucasians. 58 In 
addition, a thorough examination of other measures of 
obesity is needed to better understand the complexities 
of obesity risk across different populations. 3 Such studies 
will be enhanced by more detailed phenotypic charac- 
terisation of study cohorts that includes body fat distri- 
bution, as this is known to be a significant risk factor for 
other chronic diseases, such as cardiovascular disease 
and type II diabetes. 59 
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